Discrete element model for crack propagation in brittle materials          下载免费PDF全文

Ba Danh Le  Georg Koval  Cyrille Chazallon 《国际地质力学数值与分析法杂志》2016,40(4):583-595

We propose a discrete element model for brittle rupture. The material consists of a bidimensional set of closed‐packed particles in contact. We explore the isotropic elastic behavior of this regular structure to derive a rupture criterion compatible to continuum mechanics. We introduce a classical criterion of mixed mode crack propagation based on the value of the stress intensity factors, obtained by the analysis of two adjacent contacts near a crack tip. Hence, the toughness becomes a direct parameter of the model, without any calibration procedure. We verify the consistency of the formulation as well as its convergence by comparison with theoretical solutions of tensile cracks, a pre‐cracked beam, and an inclined crack under biaxial stress. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method          下载免费PDF全文

A. Verde  A. Ghassemi 《国际地质力学数值与分析法杂志》2016,40(6):865-886

An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading          下载免费PDF全文

K. Duan  C.Y. Kwok  M. Pierce 《国际地质力学数值与分析法杂志》2016,40(8):1150-1183

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Modeling the 2D behavior of dry‐stone retaining walls by a fully discrete element method          下载免费PDF全文

James J. Oetomo  Eric Vincens  Fabian Dedecker  Jean‐Claude Morel 《国际地质力学数值与分析法杂志》2016,40(7):1099-1120

The dry‐stone retaining walls (DSRW) have been tipped as a promising solution for sustainable development. However, before recently, their behavior is relatively obscure. In this study, discrete element method (DEM) approach was applied to simulate the plane strain failure of these walls. A commercial DEM package (PFC2D™) was used throughout this study. The authors used a fully discrete approach; thus, both the wall and the backfill were modeled as discrete elements. The methodology for obtaining the micromechanical parameters was discussed in detail; this includes the three mechanical sub‐systems of DSRWs: wall, backfill and interface. The models were loaded progressively until failure, and then the results were compared with the full‐scale experimental results where the walls were loaded, respectively, with hydrostatic load and backfill. Despite its complexity and its intensive calculation time, DEM model can then be used to validate a more simplified approach. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Modelling landslides in unsaturated slopes subjected to rainfall infiltration using material point method          下载免费PDF全文

Samila Bandara  Alessio Ferrari  Lyesse Laloui 《国际地质力学数值与分析法杂志》2016,40(9):1358-1380

This paper presents a dynamic fully coupled formulation for saturated and unsaturated soils that undergo large deformations based on material point method. Governing equations are applied to porous material while considering it as a continuum in which the pores of the solid skeleton are filled with water and air. The accuracy of the developed method is tested with available experimental and numerical results. The developed method has been applied to investigate the failure and post‐failure behaviour of rapid landslides in unsaturated slopes subjected to rainfall infiltration using two different bedrock geometries that lie below the top soil. The models show different failure and post‐failure mechanisms depending on the bedrock geometry and highlight the negative effects of continuous rain infiltrations. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Comparison of runoff characteristics of two adjacent basins in a tropical rainforest using a modified hydrologic cycle model with outflow          下载免费PDF全文

Taeko Wakahara  Katsushige Shiraki  Masakazu Suzuki 《水文研究》2014,28(3):509-520

We propose a new runoff model including an outflow process that was applied to two adjacent basins (CL, TL) located in Lambir Hills National Park in north‐central Sarawak, Malaysia. Rainfall, runoff, topography, and soil layer thickness were observed. About 19% of annual runoff was observed in the CL basin (21.97 ha), whereas about 46% was observed in the TL basin (23.25 ha). It was inferred that the CL basin has an outflow because of low base flow, small runoff peak, and excessive water loss. By incorporating the outflow process into the HYdrological CYcle MODEL, good agreement between the data generated by the model and that observed was shown, with the exception of the data from the rainless period. Then, the fitting parameters for each basin were exchanged, except for the outflow parameter, and the characteristics of each basin were compared by calculating virtual runoff. As a result, the low base flow of the CL basin was estimated by the movement of the rainwater that escaped from the basin as deep percolation or lateral flow (11% of rainfall). The potential of the CL basin for mitigating flood and drought appeared to be higher than that of the TL basin. This is consistent with the topographic characteristics of the CL basin, which has a gentler slope than the TL basin. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Numerical and analytical observations on long and infinite slopes     

D. V. Griffiths  Jinsong Huang  Giorgia F. deWolfe 《国际地质力学数值与分析法杂志》2011,35(5):569-585

Interest in the mechanics of landslides has led to renewed evaluation of the infinite slope equations, and the need for a more general framework for estimating the factor of safety of long and infinite slopes involving non‐homogeneous soil profiles. The paper describes finite element methods that demonstrate the potential for predicting failure in long slope profiles where the critical mechanism is not necessarily at the base of the soil layer. The influence of slope angle is also examined in long slopes, leading to some counter‐intuitive conclusions about the impact of slope steepness on the factor of safety. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Fracture mechanics analysis of multiple cracks in anisotropic media     

Chia‐Huei Tu  C.‐S. Chen  T.‐T. Yu 《国际地质力学数值与分析法杂志》2011,35(11):1226-1242

This paper presents a single‐domain boundary element method (BEM) for linear elastic fracture mechanics analysis in the two‐dimensional anisotropic material. In this formulation, the displacement integral equation is collocated on the un‐cracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. A special crack‐tip element was introduced to capture exactly the crack‐tip behavior. A computer program with the FORTRAN language has been developed to effectively calculate the stress intensity factors of an anisotropic material. This BEM program has been verified having a good accuracy with the previous researches. Furthermore, by analyzing the different anisotropic degree cracks in a finite plate, we found that the stress intensity factors of crack tips had apparent influence by the geometry forms of cracks and media with different anisotropic degrees. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Characterizing and classifying urban watersheds with compositional and structural attributes     

Joseph M. Delesantro  Jonathan M. Duncan  Diego Riveros-Iregui  Joanna R. Blaszczak  Emily S. Bernhardt  Dean L. Urban  Lawrence E. Band 《水文研究》2021,35(9):e14339

Current land-use classifications used to assess urbanization effects on stream water quality date back to the 1980s when limited information was available to characterize watershed attributes that mediate non-point source pollution. With high resolution remote sensing and widely used GIS tools, there has been a vast increase in the availability and precision of geospatial data of built environments. In this study, we leverage geospatial data to expand the characterization of developed landscapes and create a typology that allows us to better understand the impact of complex developed landscapes across the rural to urban gradient. We assess the ability of the developed landscape typology to reveal patterns in stream water chemistry previously undetected by traditional land-cover based classification. We examine the distribution of land-cover, infrastructure, topography and geology across 3876 National Hydrography Dataset Plus catchments in the Piedmont region of North Carolina, USA. From this dataset, we generate metrics to evaluate the abundance, density and position of landscape features relative to streams, catchment outlets and topographic wetness metrics. While impervious surfaces are a key distinguishing feature of the urban landscape, sanitary infrastructure, population density and geology are better predictors of baseflow stream water chemistry. Unsupervised clustering was used to generate a distinct developed landscape typology based on the expanded, high-resolution landscape feature information. Using stream chemistry data from 37 developed headwater catchments, we compared the baseflow water chemistry grouped by traditional land-cover based classes of urbanization (rural, low, medium and high density) to our composition and structure-based classification (a nine-class typology). The typology based on 22 metrics of developed landscape composition and structure explained over 50% of the variation in NO₃⁻-N, TDN, DOC, Cl⁻, and Br⁻ concentration, while the ISC-based classification only significantly explained 23% of the variation in TDN. These results demonstrate the importance of infrastructure, population and geology in defining developed landscapes and improving discrete classes for water management.  相似文献   

A preliminary estimate of how stream water age is influenced by changing runoff sources in the Nagqu river water shed,Qinghai-Tibet Plateau     

Yuheng Yang  Baisha Weng  Denghua Yan  Xiaoyan Gong  Yanyu Dai  Yongzhen Niu 《水文研究》2021,35(10):e14380

  相似文献